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Abstract. - Codeposition of impurities during the growth of a vicinal surface leads to an 
impurity concentration gradient on the terraces, which induces corresponding gradients in the 
mobility and the chemical potential of the adatoms. Here it is shown that the two types of 
gradients have opposing effects on the stability of the surface: Step bunching can be caused by 
impurities which either lower the adatom mobility, or increase the adatom chemical potential. 
In particular, impurities acting as random barriers (without affecting the adatom binding) 
cause step bunching, while for impurities acting as random traps the combination of the two 
effects reduces to a modification of the attachment boundary conditions at the steps. In this 
case attachment to descending steps, and thus step bunching, is favored if the impurities bind 
adatoms more weakly than the substrate. 



Introduction. Step bunching is a morphological instability of a vicinal crystal surface, in 
which a regular train of equally spaced steps separates into regions of high step density - the 
step bunches - and large flat terraces. The process can be driven energetically by an attractive 
step-step interaction, or by a variety of kinetic mechanisms, which all share the common 
feature of breaking the symmetry between the ascending (upper) and descending (lower) step 
bordering the vicinal terrace [[l] . In growth or sublimation, the symmetry breaking is provided 
by the different kinetic rates for the attachment and detachment of adatoms [0 (or some other 
species required for growth H) at the upper and the lower step; step bunching occurs under 
growth if atoms attach preferentially to the descending step. For electromigration-induced step 
bunching, the asymmetry is introduced by the electric field, and the step train is unstable if 
the adatom motion is biased in the down-step direction Q . 

It has been appreciated for a long time that in many cases step bunching must be attributed 
to the presence of impurities [|| . The traditional view is that impurities pin the steps [^) . Once 
a step is slowed down relative to its neighbors, more impurities accumulate in front of it and 
delay it even further, leading to a feedback mechanism which drives the instability 0. A 
different kind of impurity-mediated step bunching was suggested in recent work on Sii-^Cy 
layers grown on Si(100) by molecular beam epitaxy, in which C plays the role of a codeposited 
impurity Q. The key observation is that different parts of the vicinal terrace have been 
exposed to the impurity flux for different durations. Therefore the impurity concentration is 
smallest on the freshly created part near the descending step, and largest near the ascending 
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(a) (b) 

Fig. 1 - Adatom energy landscapes used in this work, (a) Random barriers: Impurities modify the 
adatom diffusion barrier by AEd while leaving binding energies unchanged, (b) Random traps: The 
adatom binding energy and the diffusion barrier for escape from an impurity site is modified by AE, 
while the transition state energies remain unchanged. 



step. To the extent that the impurities couple to the energetics and kinetics of the adatoms on 
the terrace, this causes corresponding gradients in the adatom chemical potential and mobility 
which break the symmetry between ascending and descending steps, and hence may lead to 
step bunching ( T ). 

For the SiC system, the experimentally observed step bunching could be reproduced in 
simulations in which the Si-C binding was assumed to be weaker than the Si-Si binding. This 
was interpreted in terms of an increase of the adatom mobility due to the impurities: The low 
concentration of impurities near the descending step was argued to lead to an accumulation 
of adatoms in these low-mobility regions, and hence to a preferential attachment to the de- 
scending step. However, the adatom flux onto a step depends not only on the adatom density 
gradient, but also on the adatom mobility, which is lower near the descending step. The ex- 
plicit calculations presented below show that the latter effect overcompensates the increase in 
the adatom concentration gradient. Impurities which increase the adatom mobility are found 
to stabilize the step train, while step bunching is induced if the adatoms are slowed down by 
the adsorbates. 

On the other hand, the chemical potential gradient induced by the impurities acts in 
the opposite direction. For impurities which bind the adatoms more strongly than the clean 
substrate, and which would therefore be expected to lower the adatom mobility, the chemical 
potential is decreased near the ascending step edge, where the impurities accumulate. This im- 
plies an uphill force on the adatoms, which, as is well known from studies of electromigration- 
induced step bunching, stabilizes the step train Q). Similarly, for impurities that bind more 
weakly than the substrate, the chemical potential gradient is destabilizing and the mobility 
gradient stabilizing. The net outcome of the two competing effects can be determined only if 
the modification of the adatom potential energy landscape caused by the impurities is precisely 
specified. Two limiting cases will be considered in detail: Random barriers which modify only 
the adatom mobility, and random traps for which binding energies and diffusion barriers are 
modified by the same amount (see Fig.|l|). Random barriers and random traps are standard 
models in the theory of diffusion in disordered media ]Tc[ |. 

Model and general solution. Figure [2] illustrates the geometry employed in the calculation. 
I consider a train of straight steps with spacing I. The deposition flux is F and the impurity 
flux F' . Impurities arc immobile, they do not desorb, and they are incorporated into the 
crystal when a step moves over them. The steps move with speed v — Fl. The exposure time 
at a distance x from the descending step is x/v, hence the stationary impurity coverage profile 



( 1 )A similar mechanism for step equalization was proposed in Ref. 
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Fig. 2 - Schematic of the unperturbed step train and the linear impurity concentration profile (|l|) 



is 

6{x) =F'x/v = 4>x/l, (1) 

where <j) — F' / F is the flux ratio. The spatial variation of the impurity concentration implies 
a corresponding variation of the effective chemical potential fi c s{x) and the adatom diffusion 
coefficient D(x), which will be specified later. Assuming that the adatom concentration n{x) 
adapts rapidly to changes in the step spacing and in the impurity profile, it can be computed 
from the inhomogeneous, stationary diffusion equation |]Tlf 



—D(x) — 
dx dx 



dn dfi cS 
— + fin 



dx dx 



F = 0, (2) 



where fj — l/k^T. This is supplemented by boundary conditions for the mass fluxes j_ and 
j + to the descending (x — 0) and ascending (x = I) step 0,0], 

j— = £>(0)[n'(0) + /3n(0)/4 g (Q)] = fe_n(0) (3) 

j+ = -D(l)[n'(l) + (3n(l)fi' eS (l)} = k + n{l), (4) 

The attachment rates , fc+ are chosen such that the attachment probability is symmetric in 
the absence of impurities. Specifically, I will consider two types of boundary conditions: Type 

I with fc_ = fc+ = k (attachment rates independent of the impurity concentration) and type 

II with k-/D(0) = k+/D(l) = A -1 (attachment rate proportional to the adatom mobility at 
the step). 

The mass fluxes j_ and j+ govern the dynamics of the vicinal surface. The two are related 
through mass conservation, j_ + j+ = Fl. To probe the stability of the uniform step train, 
consider a period-2 perturbation in which the length of every second terrace is increased by 
an amount e and every second terrace length is decreased by e. In the absence of impurities 
this does not affect the speed of the steps, since the attachment rates k± are symmetric, and 
hence the total flux feeding each step remains Fl. The impurity profile associated with the 
perturbed step train is therefore still given by (Q). When the coupling of the impurities to 
the adatom concentration is turned on, the larger terraces either shrink, restoring the uniform 
step train, or grow, leading to step doubling and, eventually, to step bunching. The large 
terraces shrink, if the speed of the corresponding ascending step is larger than that of the 
descending step, i.e. if 



3+(l + e) + J- (J - e) > U(l - e) + J- (J + e)- 



(5) 
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In the limit e — > this becomes dj+/dl > dj-/l or, using mass conservation, 

dj-(l)/dl<F/2, (6) 

which has to be evaluated at fixed 6{x), i.e. without taking into account the Z-dependence of 
the impurity concentration gradient. When the stability criterion (^3j) is violated, the growth 
rate of the perturbation determines the time scale for step bunching, which is given by 

r = (4dj-/dl-2F)- 1 . (7) 

Following |ll]], the mass flux j_ to the descending step is obtained from Eqs.(||,|],||) in the 
general form 

3- = Fl ( Ml t^+ . ) , (8) 
\M + l/k + + l/k-J 

M v = r u f dx x v e l3 ^ a{x) D{x)- 1 (9) 



where 



and 

fc_ = e-^^fe-, ~k+ = e-^Wfef. (10) 

Random barriers. Consider first the case where the impurities affect only the mobility of 
the adatoms. This corresponds to the random barrier energy landscape illustrated in Fig.[l] 
(a), where the binding energies remain unaffected while the diffusion barriers are modified 
by an amount AEp, which can be positive (as in Fig.[l] (a)) or negative. An exact analytic 
expression for the effective diffusion coefficient in two dimensions is not available for the 
random barrier model but some general conclusions can be drawn from Eqs.(||,^|), which 
are to be evaluated with /x e g = 0. Since D(x) is a monotonic function of x, we have that 
l/D(0) < M < l/D(l), l/2D(0) < Mi < l/2D(l) and Mi/M > 1/2 for AE D > 0, and the 
converse inequalities for AEd < 0. Using these relations it is straightforward to prove that 
attachment is primarily to the descending step (j_ > Fl/2 > j+) when AEd > 0, and to the 
ascending step when AEd < 0, for both types of boundary conditions. 

For the explicit calculation of j_ I use the expression p3,n3] 



n( \ A) Da / -i -i \ 

{X> ~ l + (ef AE ° -l)6{x) ~ 1 + { 1 

which is exact for the one- dimensional random barrier model Here Do is the diffusion 
coefficient on the clean surface, and b = (f>(e l3AED — l)/l is a dimensionless parameter describing 
the strength and the sign of the mobility gradienijpj. For type I boundary conditions the flux 
to the descending step reads 

. = Fl 1 + l/2\ a + M 2 /3A 

J ~ 2 1 + //2A + M 2 /4A ' [ ' 



where Ao = Do/k. Taking the derivative of ( |l 2|) at fixed b, one finds that the stability 
criterion (||) is satisfied (violated) when b < (b > 0). Thus step bunching is induced by 
impurities which slow down adatom diffusion (AEjj > 0). The instability is directly linked to 
the preferential feeding of the steps from above, i.e. the stability criterion (^|) is equivalent to 
j- < Fl/2. This need not be true in general. 



( 2 )Note that bl > —<f> with <f> = F' / F <C 1, hence D(x) is positive and finite everywhere. 
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Fig. 3 - Adatom density profile ( |l3| ) for i/Ao = 10 and bl = 2, -1/2 and 0, respectively. The adatom 
density lias been scaled by the overall factor F/Dq. 



The corresponding adatom density profile is given by 

n(x) = \A(X + x) + ±(Ab - l)x 2 - ^bx 3 j , (13) 

where A = j_ / F = kn(0)/F. The examples depicted in Fig|^ show how the density maximum 
shifts towards the ascending (descending) step for b > (b < 0), as would be expected intu- 
itively. Consequently the density gradient is enhanced near the ascending (descending) step. 
As was mentioned already, this effect is however overcompensated by the spatial dependence 
of the adatom mobility. The boundary values of the adatom density vary in the opposite 
direction to the density gradients, so that the mass flux is predominantly to the descending 
(ascending) step for b > (b < 0). 

For type II boundary conditions the flux to the descending step is given by 

. _ n 1 + l/2A + bl + bl 2 /3\ 

J ~~ 2 1 + l/2\ Q + bl/2 + bl 2 /4\ ' 1 ' 



which behaves similar to (12): Attachment is primarily to the descending (ascending) step and 
the step train is unstable (stable) when b > (b < 0). This shows that step bunching caused 
by random barrier impurities is a robust phenomenon which is independent of the detailed 
model assumptions. 

Random traps. In general, the influence of the impurity-induced chemical potential gradi- 
ent has to be taken into account as well. The effective chemical potential [i c g is obtained from a 
thermodynamic argument. We assume that the impurities modify the adatom binding energy 
by an amount AEb, AEf, > corresponding to stronger binding. The equilibrium adatom den- 
sity in a region with impurity concentration 8 is then given by no(0) = no(Q)(l — + 8e@ AEb ), 
since the occupation of impurity sites is enhanced or suppressed by the Boltzmann factor 
e /3AB b \yriting no(9) — no(0)e _/3AIoft and inserting the linear impurity profile (||), we obtain 



Ves(x) = -k B T dn(l + fx) 



(15) 
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where / = (e /3A - Eb — l)((f>/l) is the analogue of b in (pi]). The effective force entering (g) is 
then given by —fid^ c s/dx — f/(l + fx), which points uphill (stabilizing the step train Q) 
when AEf, > 0, and downhill (destabilizing the step train) when AEf, < 0. 

In contrast to the impurity-induced mobility gradient, the gradient in the chemical po- 
tential cannot occur in isolation (at constant D(x)), because an energy landscape in which 
the impurities modify the adatom binding energy without affecting the diffusion barriers is 
not conceivable. A simple yet realistic situation where both effects are present simultaneously 
is provided by the random trap model, illustrated in Fig.^(b). In this model it is assumed 
that the transition states between diffusion sites remain unaffected by the impurities, so that 
the binding energies and the diffusion barriers (for jumps away from an impurity site) are 
modified by the same amount, AEi, = AEq = AE. The effective diffusion coefficient is then 
given exactly by ( |ll| ) in all dimensions Jic|| . Combining (15) and ( |i"f| ) with b = /, the inte- 



grand e^ Mcff D(x)^ 1 in (j^) is seen to become constant. Hence Mo = I/Dq, Mi — 1/2Dq and 
(||) reduces to the familiar expression j| for the clean surface, but with modified attachment 
rates fc_ = fc_, fc+ = k + (l + bl). 

This is a consequence of the fact that for unbiased potential landscapes, in which the jump 
rates away from a given site are everywhere symmetric, the inhomogeneous diffusion equation 
(|2|) can be written as [p|,|l5|] 

^[D(x)n(x)]+F = 0, (16) 

which implies that h(x) = D(x)n(x) satisfies a diffusion equation with constant coefficients 
and standard boundary conditions D(0)n'(0) = fc_n(0), D(l)n'(l) = — k + n(l). For type 
II boundary conditions k + /D(l) = fc_/_D(0) so that attachment remains symmetric in the 
presence of impurities, and the two competing impurity effects precisely cancel. In contrast, 
for type I boundary conditions the impurities are seen to induce a preference for attachment 
at the ascending step for AE > 0, and at the descending step for AE < 0. This implies 
a tendency towards step bunching for AE < 0. The effect is however quite feeble, since 
\bl\ < 1 for AE < 0. 

Summary. In conclusion, I have described a novel mechanism through which codeposited 
adsorbates may destabilize a growing vicinal surface. The kinetic and energetic couplings 
between adsorbate atoms and adatoms were shown to have competing effects. Impurities 
which slow down the adatom diffusion without affecting adatom binding energies (random 
barriers) generically cause step bunching. When the impurities act as random traps, and 
provided the attachment rates at the steps are not modified by the impurities (type I boundary 
conditions) , the net result of the two effects was found to be destabilizing for impurities that 
bind adatoms more weakly than the substrate. This is consistent with the simulations of SiC 
growth |8| that inspired the present study. However, the precise cancellation of the two effects 
for type II boundary conditions also suggests that predicting the stability of the surface may 
be difficult if the details of the adsorbate-adatom-interaction are not known. 

The observation that barrier-like and trap-like impurities may have qualitatively different 
effects on the stability of a growing surface, because they affect the symmetry of the surface 
diffusion process in different ways, was made previously in the context of growth on singular 
surfaces [jj"6[ . It is interesting to note that also in this case the barrier-like impurities induce 
a downhill diffusion bias, favoring attachment to the descending step, through a mechanism 
that however does not involve any impurity concentration gradient. 

Future work should address the competition between the impurity-induced instability and 
the stabilizing effect of conventional step edge barriers Going beyond the linear stability 
analysis presented here will be difficult because the dynamics becomes nonlocal in time when 
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the impurity profile is nonstationary j7|. Further underpinning for the proposed mechanism 
from KMC simulations would therefore be highly desirable. 

* * * 
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